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摘 要 土壤 养分 影响 着 土壤 的 质量 ， 也 影响 着 植被 、 农 作物 等 的 生长 。 为 快速 准确 地 估 测 艾 比 湖 流域 土壤 
养分 状况 ,选择 艾 比 湖 流域 精 河 县 作为 研究 区 ， 以 精 河 县 内 不 同 地 表 和 覆盖 类 型 土壤 为 研究 对 象 ， 基 于 实地 采 
集 的 75 个 土壤 样品 的 室内 ASD FieldSpec3 实测 光谱 数据 和 3 种 光谱 变换 形式 , 利用 10 nm 间隔 重 采 样 进行 去 
噪 处 理 ， 再 结合 多 元 逐步 回归 法 (SMLR)、 偏 最 小 二 乘法 回归 法 (PLSR)、 人 工 神经 网 络 法 (ANN) 分 别 建立 土壤 
养分 预测 模型 ， 以 探索 最 优 模 型 .结果 表明 : 土壤 实测 光谱 的 一 阶 微分 、 二 阶 微分 变换 形式 能 显著 提高 光谱 与 
土壤 养分 之 间 的 相关 性 ， 尤 其 是 一 阶 微分 变换 与 土壤 有 机 质 和 全 和 氮 的 相关 性 最 高 分 别 达 0.87 和 0.91， 光 谱 变 
换 技 术 能 显著 增强 土壤 养分 与 高 光谱 之 间 的 敏感 度 ， 达 到 更 好 的 建 模 效果 ; SMLR、PLSR 和 ANN 这 3 种 模型 
都 具有 良好 的 预测 能 力 ， 其 中 , ANN 建立 的 模型 预测 效果 最 好 ,二 阶 微 分 变换 的 ANN 模型 对 有 机 质 、 全 和 氮 的 
硕 测 决 定 系数 (R”) 分 别 为 0.886 和 0.984, 均 方 根 误差 (RMSE) 分 别 为 2.614 和 0.147, PLSR 次 之 ; 全 氮 的 预测 效 
明显 优 于 有 机 质 的 预测 效果 ,说 明 高 光谱 和 全 和 氮 之 间 的 敏感 性 更 高 . 总 体 来 说 ,光谱 二 阶 微 分 变换 形式 的 人 
工 神经 网 络 模型 可 以 最 精确 稳定 地 完成 土壤 养分 含量 的 快速 预测 ， 能 够 实现 艾 比 湖 流域 的 土壤 养分 空间 分 布 
状况 和 动态 变化 特征 的 动态 监测 。 

关键 词 ”土壤 养分 高 光谱 遥感 ”多 元 逐步 回归 偏 最 小 二 乘法 “人工 神经 网 络 艾 比 湖 流 域 
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Abstract The soil nutrient affects soil quality, vegetation type, crops growth and yield. To rapidly and accurately determine 
Soil nutrient contents, an indoor spectral data (measured by ASD FieldSpec3) of 75 soil samples of Jinghe County of Ebinur 
Lake Basin were analyzed. Then the collected data were processed at resampling interval of 10 nm to suppress noise. Soil 
nutrient hyper-spectral forecast models were used to forecasts soil nutrient contents in three transformation conditions. The 
performance of the models was evaluated based on stepwise multiple linear regression (SMLR) analysis, partial least squares 
regression (PLSR) analysis and artificial neural network (ANN) analysis and the optimal model determined by comparison. 
The results showed that the transformation of first-order and second-order differential dramatically enhanced correlation 


between spectroscopy data and soil nutrient content. Specifically, the first-order differential of soil spectroscopy had a good 
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correction with soil nutrient content. The correlation coefficients for organic matter and total nitrogen were 0.87 and 0.91, 


respectively. In conclusion, spectral transformation technique increased the sensitivity of high Spectral data to soil nutrient 


change, and it produced far better forecasting results. Although all the three models had good predictive ability, ANN model 
had the best predictive effect, followed by PLSR model. The ANN model estimation test based on the second-order differential 
of Spectroscopy data with independent datasets from different soil samples respectively produced R2 and RMSE values of 


0.885 and 0.984 for organic matter and 2.614 and 0.147 for total nitrogen. The prediction effect of total nitrogen was obviously 


better than that of organic matter. This indicated that the sensitivity of soil hyper-spectral reflectance to the soil total nitrogen 


content was much better. Overall, the ANN model based on the second-order differential of Spectroscopy data rapidly and 


precisely predicted soil nutrients contents. It was beneficial for monitoring Spatial distributions and dynamic changes of soil 


nutrients in Ebinur Lake Basin. 
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土壤 是 陆地 图 层 能 生长 植物 的 下 松 表层 ,能 
持续 、 全 面 地 为 植物 生长 提供 所 需 的 水 、 肥 、 气 、 热 
等 基本 要 素 …。 土 壤 养 分 是 土壤 系统 的 重要 组 成 ， 影 
响 着 土壤 的 理化 性 质 ， 是 制约 植被 生长 和 土地 生产 
能 力 的 关键 因素 门 。 土 壤 有 机 质 和 全 氮 合 量 是 植物 生 
长 的 主要 养分 , 同时 也 是 土壤 微生物 生命 活动 的 主 
要 能 量 来 源 " 1 其 含量 能 直接 体现 土地 的 生产 能 力 。 
因此 土壤 有 机 质 、 全 氮 含 量 及 其 动态 平衡 是 反映 土壤 
健康 和 土壤 质量 的 重要 指标 。 土 壤 养 分 含量 的 快速 
精准 测定 是 现代 精细 农业 发 展 的 重要 前 提 户 。 传 统 土 
二 养分 含量 的 测定 往往 通过 试验 分 析 获 得 ,试验 结 
果 相 对 精准 , 但 试验 分 析 所 需要 的 经 济 成 本 和 时 间 
成 本 较 高 ,不 适宜 大 尺度 范围 土壤 养分 的 测定 站。 

土壤 反射 光谱 与 土壤 理化 性 质 之 间 的 密切 关系 ， 
为 土壤 组 分 和 属性 的 穆 感 研究 提供 了 可 能 “"。 大 量 
的 研究 围绕 土壤 反射 光谱 与 土壤 有 机 质 、 气 素 、 水 
分 、 重 金属 等 进行 ， 并 取得 了 不 错 的 反 演 效果 ”)。 
国内 外 很 多 学 者 都 使 用 近 红外 光谱 对 土壤 组 分 反 演 
进行 建 模 。Hummel 等 "分 析 了 近 红 外 光谱 与 土壤 
有 机 质 之 间 的 关系 ,利用 多 元 逐步 线性 回归 建立 有 
机 质 预测 模型 。 张娟娟 等 "建立 了 基于 不 同 光 谱 指 
数 的 土壤 有 机 质 光 谱 估算 模型 。 任 红 艳 等 "综合 利 
用 可 见 光 - 近 红外 光谱 、 近 红外 /中 红外 光谱 结合 含 
最 小 二 乘 回归 分 析 成 功 地 预测 了 不 同 区 域 的 土壤 有 
机 左 、 全 氮 的 合 量 。 李 伟 等 (同时 利用 偏 最 小 二 乘 
法 和 人 工 神经 网 络 法 建立 了 黑土 的 碱 解 氮 、 速 效 磷 、 
速效 钾 含 量 的 预测 模型 ， 发 现 人 工 神经 网 络 法 有 一 
个 更 好 的 建 模 效 果 。 

高 光谱 遥感 凭借 其 高 分 辩 率 等 优点 ， 能 够 更 好 
地 反 演 土壤 组 分 特征 ,更 有 利于 土壤 养分 的 遥感 估 
算 。 关 于 高 光谱 遥感 监测 土壤 的 研究 , 大 多 借鉴 红 
外 光谱 遥感 技术 研究 的 经 验 , 在 原 有 的 遥感 光谱 模 
型 基础 上 进行 改进 。 刘 华 等 (小 王 维 等 "4、 贺 军 亮 
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等 ("都 对 不 同 环境 背景 下 的 土壤 重金 属 高 光谱 特征 
进行 了 较为 细致 的 研究 ， 均 得 到 不 错 的 效果 。 陈 东 
强 等 六 基于 高 光谱 技术 研究 了 准噶尔 盆地 人 工 林 全 
盐 含量 ， 进 而 了 解 土地 退化 程度 。 李 民 赞 (对 不 同 
环境 背景 下 的 土壤 有 机 质 高 光谱 特征 进行 了 细致 的 
研究 。 目 前 高 光谱 遥感 监测 土壤 组 成 的 常用 统计 分 
析 方 法 有 多 元 逐步 回归 分 析 握 }、 仿 最 小 二 乘 回归 分 
析 号 、 人 工 神经 网 络 分 析 扩 1;。 大 量 研究 证 明 这 3 种 
方法 在 高 光谱 遥感 建 模 上 具有 适用 性 。 新 疆 因 其 
独特 的 气候 特征 和 地 理 环 境 条 件 ， 形 成 了 典型 的 
山地 -绿洲 -荒漠 系统 中 。 对 于 这 种 特殊 景观 的 土 
坏 属 性 高 光谱 研究 并 不 多 ， 且 多 集中 于 土壤 盐 渍 化 
层面 ("而 关于 新 疆土 壤 养 分 的 高 光谱 研究 尚 不 
多 见 。 使 用 目前 已 有 的 高 光谱 经 验 模 型 无 法 精确 反 演 
新 疆 地 区 的 土壤 组 分 和 属性 状况 外。 本 研究 以 新 疆 
艾 比 湖 流域 不 同 地 表 覆 盖 类 型 的 土壤 养分 为 研究 对 
象 , 通过 光谱 变换 技术 进行 优化 处 理 ， 选 取 多 元 逐步 
回归 分 析 (stepwise multiple linear regression, SMLR)、 
偏 最 小 二 乘法 分 析 (partial least squares regression, 
PLSR) 和 人 工 神经 网 络 分 析 (artificial neural network， 
ANN)3 种 方法 进行 高 光谱 模型 建立 ， 试 图 通过 分 析 
3 种 模型 的 稳定 性 和 准确 性 ， 探 求 艾 比 湖 流域 土壤 
养分 高 光谱 反 演 的 适宜 性 和 可 能 性 ， 以 期 得 到 能 
演 土壤 养分 空间 分 布 状况 和 动态 变化 特征 的 最 佳 模 
型 , 为 艾 比 湖 流域 人 工 绿 洲 的 稳定 、 农 业 发 展 以 及 
荒漠 化 治理 提供 科学 依据 。 
1 材料 与 方法 
1.1 土壤 样品 采集 

本 研究 以 天 山北 坡 艾 比 湖 流域 的 精 河 县 
(81?"46'83"51'E, 44°02~45°10'N) 为 试验 点 。 该 区 属于 上 典 
型 的 干旱 荒漠 型 大 陆 性 气候 ,年 均 降 水 仅 有 90.9 mm 
而 年 均 蒸发 量 高 达 1 423.9 mme 年 均 气 温 为 8 CC， 
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多 大 风 天 气 ， 是 西北 沙尘暴 重要 源 区 之 一 。 精 河 县 
地 形 复 杂 ， 三 面 环 山 (婆罗 科 努 山 、 科 古琴 山 、 阿 拉 
套 山 )， 东 北面 与 准噶尔 盆地 相连 接 , 形成 了 由 南 和 
西南 向 北 倾斜 的 艾 比 湖 为 底 中 心 的 盆地 , 平原 多 为 
山 前 坡 积 洪 积 平原 ,广阔 起 伏 ， 上 面 广 泛 分 布 着 砾 
质 戈壁 、 风 积 沙丘 和 盐 碱 沼泽 ， 是 典型 的 山地 -绿洲 
- 苇 漠 系统 。 土 壕 类 型 众多 ,地表 覆盖 类 型 复杂， 包 
括 厌 棕 漠 土 、 灰 漠 土 、 草 旬 土 、 沼 泽 士 、 潮 土 、 灌 
淤 士 、 风 沙土 和 盐 士 8 种 土壤 类 型 。 于 2015 年 8 月 ， 
选取 梭 梭 (Haloxylon ammrodenradromJ、 芒 革 (Pjragmaites 
australis)、 胡杨 (Populus euphratica)、 构 杞 (Lyciumbarbarum)、 
棉花 (4nemone vitifolia)、 沙 地 又 (Morus alba)、 罗 布 
麻 (4pocynum venetum) 等 15 种 不 同 地 表 覆 盖 类 型 ， 


样 地 基本 信息 见 表 1。 样 地 选取 不 仅 考虑 具有 代表 
性 的 覆盖 类 型 ， 还 考虑 了 空间 分 布 。 本 研究 的 样 地 
选取 位 于 精 河 县 人 工 绿 洲 规划 区 和 生态 修复 区 , 集 
中 于 精 河 县 中 部 和 东北 部 平原 地 区 ， 围绕 G312 干 
线 ,避免 过 多 的 人 为 扰动 ,大 多 选取 在 距离 公路 起 
过 1 km 处 ， 见 图 1。 在 各 样 地 进行 土壤 剖面 调查 ,每 
个 样 地 控 取 5 个 100 cm 深 的 土壤 剖面 (重复 试验 )， 
使 用 100 cm 环 刀 分 5 层 分 层 取 样 ， 再 用 小 铲子 每 个 
样 取 土 1 kg， 同 时 使 用 手持 式 GPS 仪 对 每 个 采样 点 
的 经 纬度 和 海拔 进行 记录 。 土 样 经 自然 风干 完全 后 ， 
研磨 、 去 除 砂 粒 、 石 块 和 植物 残 体 , 用 0.25 mm 孔 
际 的 筛子 进行 筛选 。 将 处 理 好 的 土 样 分 为 3 份 , 分 
别 用 于 光谱 分 析 、 化 学 测量 和 保存 备用 。 


表 1 
Table 1 


试验 调查 样 地 的 基本 信息 


Information of the investigated sample plots 


样 地 编号 


经 纬度 


海拔 


- 。 土壤 类 型 土地 利用 类 型 地 表 覆 盖 类 型 (主要 植物 类 型 ) 
Sample plot, Latitude and, altriude Agrotype Land-use type Land cover type (major plants) 
number longitude (m) EOP yP ° yP Pn 
J、 模 楼、 和 穴 条 、 麻 黄 、 沙 近 志 
1# 83?19 下 ， 3 盐 土 草地 | 白 刺 a 麻黄 Be 到 
44°34N a jesture and itraria tangutorum, Haloxylon ammodendron, Caragana 
korshinskii, Ephedra sinica, Calligonum arborescens 
蔡 、 赴 站 二 
83°09E, 1] 沼泽 土 水 域 本 i eli 
44°37'N Bog soil Water area hragmites australis, Halocnemum strobilaceum, Kalidium 
foliatum 
R j、 、 、 
83°29'E, 灰 漠 土 未 利用 土地 和 
3# o 369 ; Nitraria tangutorum, Haloxylon ammodendron, Tamarix 
44°33'N Gray desert soil Unused land ， 
hohenackeri, Populus euphratica 
83°41 玉 潮 土 草地 区 | 
4# > 330 罗布 床 、 上 4pocynum venetum, Phragmites australis 
44°36'N Fluvo-aquic soil Pasture land ee S 
83"34'， 风沙 土 未 利用 土地 无 植被 
5# 294 , 
44°40'N Aeolian sandy soil Unused land Unvegetated 
、 沙 拐 京 、 刺 沙 鞍 、 盐 
6# 83°12 下 ， 3 灰 棕 漠 土 耕地 有 Cn ee 盐 生 草 po 
44°37'N rl Cultivated land aloxylon ammodendron, Elaeagnus angustifolia, Salsola 
ruthenica, Halogeton glomeratus 
82°56'E, 灰 漠 土 耕地 枸杞 
7# > 324 , 
44°36'N Gray desert soil Cultivated land Lycium barbarum 
8 82°56'E, a 灰 漠 土 林地 梭 梭 
44°33'N Gray desert soil Forest land Haloxylon ammodendron 
时 沙 束 
0 82°57'E， 号 灌 淤 土 林地 、 树 柳 
44°35'N 2 Forestland opulus simonii, UlImus pumila, Elaeagnus angustifolia, 
Tamarix hohenackeri 
0# 82°57' 祷 ， a 草 旬 土 草地 沙 束 、 沙 拐 刺 
44°35N Meadow soil Pasture land Elaeagnus angustifolia, Calligonum arborescens 
这 82°45'E, a 灌 淤 土 耕地 沙 地 桑 
44°31'N Anthropogenic-alluvial soil Cultivated land Moru salba 
站 3 二 
82o44， 风沙 二 未 利用 土地 人 
2# o 372 Haloxylon ammodendron, Halocnemum strobilaceum, 
44°31N Aeolian sandy Soil Unused land i 
Kalidium foliatum 
3 82°46'E, 3 灌 淤 土 耕地 棉花 
44°31'N Anthropogenic-alluvial soil Cultivated land Anemone vitifolia 
82°45'E, 盐 土 林地 小 叶 杨 
4# 330 
44°31'N Solonchak Forest land Populus simonii 
二 82"38'， 5 灌 淤 土 耕地 枸杞 
44°33'N Anthropogenic-alluvial soil Cultivated land Lycium barbarum 


1.2 土壤 养分 合 量 测定 

将 各 取样 点 的 土壤 样品 带 回 实验 室 ， 使 用 恒温 
105 Y 干 燥 箱 烘 干 至 恒 重 ， 测 得 容重 。 将 已 去 除 杂 
质 、 烘 干 的 土壤 溶 于 蒸馏 水 中 , 用 pH 计 测 土壤 pH。 


土壤 全 赴 含 量 的 测定 使 用 残渣 - 烘 干 质量 法 1， 土壤 
有 机 质 含量 的 测定 采用 重 铬 酸 钾 容 重 - 外 加 执法 进行 
测定 9 全 氮 含 量 则 选用 最 常用 的 凯 氏 定 氮 法 MA。 表 
2 为 不 同样 点 土壤 样本 的 基本 性 质 ， 包括 土壤 容重 、 


http:/www.ecoagri.ac.cn 


1558 中 国生 态 农业 学 报 2016 


2 


45°N 


0 35 70km 人 
上 一 一- 


绕 
区 这 
Rey ROCA Ci b .党 党 


和、 这 “aX 党 
Ww 
了 | “采样 点 
2| Sampling plot 
82°E 83°E 
1 研究 区 及 采样 点 示意 图 
Fig.1 Location of the research area and distribution of 


sample plots 


pH、 总 盐 等 基本 参数 。 
1.3 ”土壤 光谱 采集 及 预 处 理 
光谱 分 析 采 用 美国 ASD 公司 生产 的 FieldSpec3 


地 物 光谱 仪 。 光 谱 仪 重 采样 间隔 为 1 nm， 共 输出 
2 151 个 波段 。 试 验 采 用 暗室 试验 ,试验 时 将 土壤 均 
匀 放 置 在 盛 样 血 内 ， 保 持 土 层 厚度 约 为 2 cm。 为 减 
少 试验 土 样 每 个 方向 的 异 质 性 ,测量 时 分 别 在 盛 样 
轩 0"、90"、180" 和 270* 方 位 各 采集 3 条 样本 线 ， 每 
个 土 样 都 采集 12 条 样本 光谱 数据 ， 进 行 算术 平均 所 
得 值 构成 曲线 则 为 光谱 曲线 , 同时 在 每 次 测试 之 前 
光谱 仪 都 需 进行 白板 标定 Ci。 光谱 采集 试验 选择 室 
内 进行 而 不 采用 野外 实地 测量 ， 是 因为 实测 光谱 数 
据 容易 受 外 界 影 响 (植被 覆盖 、 土 壤 水 分 、 大 气 条 
件 等 )， 尤其 是 研究 地 精 河 县 大 风 天 气 频 发 ,扬尘 
较 多 ,选择 室内 进行 试验 可 控 性 强 ， 易 获得 稳定 的 
光谱 数据 。 


表 2 研究 区 各 采样 点 土壤 基本 性 质 


Table 2 Basic properties of soil samples of soil sample plots 


样 地 编号 容重 总 盐 有 机 质 全 氮 
Sample plot number Bulk density (g-cm) ph Total salinity (g-kg !) Organic matter (gkg !) Totalnitrogen (gkg !) 
1# 1.160+0.028 8.090+0.355 2.40+0.324 4.468+0.114 0.191+0.090 
2# 1.519+0.031 7.730+0.075 14.78+0.870 4.252+0.072 0.114+0.012 
3# 1.196+0.070 7.620+0.219 19.70+0.957 4.980+0.077 0.420+0.015 
4# 1.375 十 0.024 8.352+0.089 27.60+0.843 3.359+0.153 0.190+0.011 
5# 1.543+0.025 7.534+0.051 3.52+0.642 0.402+0.068 0.190+0.021 
6# 1.621+0.077 7.552+0.089 4.78+0.259 0.843 土 0.037 0.038+0.008 
7# 1.752+0.070 7.466+0.078 7.10+0.442 2.991+0.137 0.114+0.016 
8# 1.604+0.050 8.202+0.390 0.88+0.148 2.294+0.128 0.114+0.013 
9# 1.157+0.065 7.958+0.165 8.10+0.447 2.244+0.095 0.191+0.015 
0# 1.156+0.056 7.744+0.151 17.30+0.949 3.827+0.168 0.114+0.010 
1# 1.407+0.068 8.428+0.328 3.08+0.228 3.763+0.098 0.190+0.014 
2# 1.041+0.044 8.718 土 0.245 2.204+0.235 1.683 二 0.093 0.114+0.019 
3# 1.442+0.067 8.522+0.286 3.80+0.265 4.385+0.186 0.152+0.004 
4# 1.805+0.039 7.796+0.154 10.34+0.820 0.457+0.025 0.114+0.010 
5# 1.593+0.049 8.542+0.067 1.88+0.228 12.200+0.298 0.649+0.028 


尽管 在 样本 光谱 数据 采集 之 前 , 土壤 已 经 过 风 
干 、 研 磨 、 过 筛 等 处 理 , 但 由 于 不 同 土壤 类 型 矿物 
组 成 等 属性 不 同 导致 土壤 样品 的 颗粒 大 小 有 细微 差 
别 ， 影响 吸光 度 的 大 小 而 干扰 光谱 信息 的 分 析 和 提 
取 Wl。 为 减少 光谱 测定 过 程 中 的 误差 , 往往 对 土壤 
原始 样本 进行 去 噪 、 重 采样 处 理 。 有 研究 表明 ， 光 
谱 不 同形 式 变 化 能 降低 背景 、 样 品 粗糙 度 等 因子 的 
干扰 FJ。 本 研究 采用 10 nm 间隔 重 采样 去 噪 ， 重 采 
样 后 能 获取 215 个 波段 数组 成 的 光谱 曲线 ， 在 更 加 
平滑 的 同时 仍 能 维持 原 光 谱 的 形状 特征 。 为 扩大 样本 
之 间 的 光谱 特征 差异 ， 更 有 效 分 析 高 光谱 特征 与 土 
壤 有 机 质 、 全 和 氮 之 间 的 关系 ， 本 研究 对 光谱 反射 率 进 
行 一 阶 微分 、 二 阶 微分 变换 ， 有 研究 表明 光谱 变换 技 


术 能 显著 提高 高 光谱 反 演 的 稳定 性 和 精确 性 请 ]。 
通过 对 所 测 土壤 光谱 反射 率 及 其 相关 变换 形式 
与 土壤 养分 含量 的 逐 波段 相关 分 析 ， 对 光谱 数据 进 
行 降 维 处 理 ， 除 去 无 效 信息 。 选 取 相 关系 数 较 大 波 
段 ， 作 为 待 建 模 因 变量 ， 避 免 了 模型 过 于 复杂 而 不 
能 快速 达到 预测 效果 。 
1.4 ”高 光谱 模型 建立 和 精度 验证 
1.4.1 ”模型 建立 
多 元 逐步 回归 分 析 (SMLR) 根 据 方差 分 析 结 果 
选取 效果 较 好 的 自 变 量 进入 方程 进行 建 模 , 能 够 将 
F 值 较 小 的 变量 剔除 ， 达 到 较 好 的 建 模 效果 ， 是 目 
前 解决 多 元 线性 问题 最 常用 的 方法 中。 本 研究 的 
SMLR 模型 是 采用 SPSS 19 处 理 得 到 的 。 偏 最 小 二 
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乘法 分 析 (PLSR) 模 型 是 目前 高 光谱 反 演 建 模 选用 较 。 验证 , 很 好 地 反映 模型 自身 的 稳定 性 ，R* 较 大 时 ， 


多 、 效 果 较 好 的 回归 模型 ， 汇 集 典 型 相关 、 主 成 分 
和 多 元 线性 回归 3 种 分 析 方 法 优点 上 2 ， 能 够 最 大 
程度 地 提取 数据 变异 信息 ， 特 别 适 合 变量 高 线性 
相关 的 数据 分 析 ; 所 建 模型 能 简化 方程 ， 并且 保 留 
大 量 信息 ， 提 高 运算 速度 P "1。 本 研究 的 PLSR 模型 
建立 是 采用 Minitab 17 拟 合 实现 的 。BP 神经 网 络 是 
人 工 神 经 网 络 (ANN) 中 运用 最 广泛 的 算法 ， 大量 研 
究 证 明 ANN 模型 中 的 BP 神经 网 络 算法 能 够 很 好 地 
利用 光谱 特征 反 演 土壤 组 分 扬 ”。ANN 模型 也 需要 
相关 系数 法 先进 行 波段 筛选 , 过 多 的 波段 会 使 网 络 
结构 复杂 ， 导 致 计算 量 大 、 收 敛 慢 等 问题 ”YY。ANN 
模型 涉及 多 套 方案 的 比较 和 簿 选 县 有 同时 处 理 线 
性 和 非 线性 关系 的 能 力 CJ。 本 研究 的 ANN 模型 建 
立 使 用 DPS 7.05。 
1.4.2 ”模型 精度 验证 

本 研究 主要 采用 决定 系数 (R*) 和 均 方 根 误 差 
(RMSE) 来 进行 模型 效果 验证 ， 回 归 模 型 的 评价 采 
用 方差 分 析 法 。R? 对 自 变量 和 因 变 量 之 间 关 系 进行 

- .- :=- 反射 率 Original spectroscopy 


a ”一 一 一 一 阶 微分 Firstrorder differential spectroscopy 
1,0 FE vires 二 阶 微分 Second-order differential spectroscopy 


相关 系数 Correlation coefficient 


模型 较 稳定 
的 程度 ， 是 反映 模型 监测 角 
RMSE 较 小 时 , 模型 精度 较 高 ， 预测 能 力 较 强 。 所 
以 模型 的 R* 越 大 , RMSE 越 小 ， 模 型 的 稳定 性 越 好 ， 
预测 精度 也 越 高 中 。 本 研究 的 R 和 RMSE 计算 以 
及 模型 拟 合 分 布 图 像 的 绘制 均 由 Microsoft Excel 
2007 实现 。 


2 结果 与 分 析 
2.1 土壤 养分 空间 分 析 及 与 光谱 的 相关 性 分 析 
由 表 2 土 样 基本 性 质 可 知 不 同 地 表 履 盖 类 型 土壤 
养分 的 含量 差异 较 大 。 2 种 养分 的 最 大 值 均 出 现在 15# 
样 地 ,土壤 有 机 质 和 全 氮 含 量 分 别 达 12.200 gkg”! 
和 0.649 gkg ; 其 次 是 3# 样 地 , 而 5#、6# 和 14# 样 
地 2 种 土壤 养分 含量 都 相对 较 低 。 
经 过 相关 系数 法 处 理 得 到 ,土壤 有 机 质 和 全 和 氮 
量 与 光谱 原始 反射 率 、 一 阶 微 分 和 二 阶 微分 之 间 
的 相关 系数 如 图 2 所 示 。 


定 ; RMSE 能 够 表达 拟 合 数据 偏离 真实 值 
E 力 的 重要 指标 ， 当 


3 


We Hh | 出 


一 1.0 
350 600 850 1100 1350 1600 1850 2100 2350 2600 350 600 850 1100 1350 1600 1850 2100 2350 2600 
波长 Wavelength (nm) 


2 土壤 有 机 质 (a) 和 全 和 所 (b) 含 


量 与 光谱 反射 率 及 其 一 阶 微分 和 二 阶 微分 的 相关 性 


Fig. 2 Correlation between soil organic matter content (a), total nitrogen content (b) and original, first-order differential and 
second-order differential soil Spectroscopy 


两 种 养分 和 原 光 谱 反 射 率 的 负 相 关 都 集中 在 
350~470 nm 处 ,并 都 在 1 880~1 890 nm、1 980~2 000 nm 
和 2 140~2 230 nm 处 有 较 好 的 正 相 关 性 , 全 氮 的 相 
关 性 略 高 于 有 机 质 ; 有 机 质 和 全 和 氮 与 一 阶 微分 光谱 
变换 形式 的 相关 性 较 高 波段 集中 于 650~840 nm 和 
880~920 nm 处 , 且 两 者 最 好 相关 性 波段 都 为 890 nm， 
有 机 质 的 相关 性 达 0.87, 全 氮 的 则 高 达 0.91; 两 种 物 
质 与 光谱 二 阶 变 换 形 式 相 关 性 较 高 波段 相对 分 散 ， 
其 中 有 机 质 在 1 410 nm 处 正 相关 最 高 达 0.81, 负 相 
关 最 高 只 有 0.70 出 现在 1 120 nm 处 。 而 全 和 氮 则 在 


770 nm、1 010 nm 和 2 430 nm 处 呈现 较 好 的 正 相关 ， 
相关 性 最 高 达 0.78, 负 相 关 较 好 处 在 740 nm、1 030 nm 
处 , 相关 性 0.79。 

从 图 2 可 明显 看 出 ,两 种 土壤 养分 与 光谱 数据 
相关 性 基本 一 致 ， 两 种 养分 实测 数据 间 的 相关 性 也 
达 0.85。 有 研究 表明 , 土壤 有 机 质 含量 与 全 氮 含 量 
之 间 存 在 密切 的 关系 ,水稻 土 全 氮 含量 可 以 用 有 机 
质 含量 乘 以 转换 系数 0.065 获取 , 干旱 地 区 的 转换 
系数 为 0.05~0.0655。 两 种 养分 相关 性 略 有 差异 是 由 
于 高 光谱 反 演 基于 各 种 不 同化 学 键 位 的 吸收 ， 而 不 
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同 物 质 组 成 导致 两 种 养分 含量 在 敏感 波段 上 出 现 一 
定 的 差异 性 0?1。 
2.2 ”多 元 逐步 回归 建 模 

对 土壤 养分 与 反射 光谱 及 其 变换 形式 之 间 进 行 
相关 分 析 , 不 同 变 换 形式 光谱 与 土壤 养分 的 较 大 相 
关系 数 所 对 应 波段 通常 作为 预测 土壤 养分 的 敏感 波 
段 ， 故 选取 各 形式 相关 系数 较 大 的 波段 作为 自 变 量 ， 
以 有 机 质 含 量 和 全 和 氮 含 量 作为 因 变 量 进行 多 元 逐步 
线性 回归 分 析 ， 获 得 光谱 原始 反射 率 、 光 谱 一 阶 微 
分 及 二 阶 微分 的 SMLR 模型 以 及 SMLR 模型 所 得 预 
测 值 与 养分 实测 值 之 间 关 系 如 表 3 所 示 。 光 谱 反 射 


率 的 一 阶 微 分 和 二 阶 微分 变换 形式 所 得 到 的 R 和 
RMSE 均 优 于 光谱 原 反射 率 的 值 ， 其 中 二 阶 微分 略 
优 于 一 阶 微分 。 有 机 质 和 全 所 二 阶 微分 所 建 模型 的 
R? 分别 达 0.803 和 0.852, RMSE 只 有 1.361 和 0.066， 
能 够 很 好 地 预测 土壤 有 机 质 和 全 氮 两 种 养分 状况 。 
一 阶 微分 建 模 所 得 及 也 能 达到 0.768 和 0.835， 对 应 
的 RMSE 也 只 有 1.42 和 0.064,， 也 能 相对 较 好 地 预 
测 土壤 养分 。 而 光谱 原 反 射 率 所 建 模型 则 相对 较 差 ， 
不 能 得 到 稳定 精确 的 预测 值 。 比 较 两 种 养分 的 预测 
情况 ,全 所 的 拟 合 效果 优 于 有 机 质 ， 高 光谱 对 全 氛 
的 敏感 性 强 于 有 机 质 。 


表 3 土壤 养分 多 元 逐步 回归 分 析 结 果 


Table 3 Test results of stepwise multiple linear regression (SMLR) model for soil nutrient contents 


多 元 逐步 回归 模型 


决定 系数 


土壤 养分 光谱 变换 形式 ee: 总 均 方 根 误 差 
2 Stepwise multiple linear Coefficient of 
Soil nutrient Transform form of spectrum Root mean Square error 
regression model determination 
NVR - 
有 机 质 、. 光谱 反射 率 7=24.193X18g0—7.679 0.486 2.113 
Organic matter Original soil spectroscopy 

NER - 一 8 人 和 

光谱 区 射 率 一 阶 微分 7=28 903.701X8oo+3.072 0.768 1.420 
First-order differential soil spectroscopy 
光谱 反射 率 二 阶 微分 =-3 554.960Xiooo+ 1 Te 
Second-order differential soil spectroscopy 74 349.124X1410-0.577 : 
NE 寸 变 
全 所 a 光谱 反射 率 7=-0.151—1.053X390t+1.517X1900 0.627 0.093 
Total nitrogen Original soil spectroscopy 

NR 村 变 二 阶 向 分 

光谱 区 射 率 一 阶 微分 二 0.170+1 611.122X%o0 0.835 0.064 
First-order differential soil spectroscopy 
光谱 反射 率 二 阶 微分 0.274- 于 

光谱 反射 率 二 阶 微分 7=0.274-1 492.812X1030 de dee 


Second-order differential soil spectroscopy 


5 786.504X610t+10 819.675X770 


XN 表示 N~N+9 (nm) 波 段 的 算术 平均 数 。X indicates the arithmetic mean of N to N+9 (nm) waveband. 


2.3 ” 偏 最 小 二 乘法 建 模 

基于 PLSR 模型 土壤 养分 实测 值 与 预测 值 的 拟 
合 分 布 如 图 3。 选 取 建 模 波段 为 自 变 量 ， 以 土壤 养分 
为 因 变 量 进行 PLSR 分 析 ， 针 对 光谱 原始 反射 率 、 一 
阶 和 二 阶 微分 变换 形式 进行 建 模 。 反 射 率 一 阶 微分 
和 二 阶 微分 变换 形式 所 建 模型 优 于 原 反射 光谱 所 建 
模型 ， 说 明光 谱 变换 技术 能 显著 提高 模型 的 预测 精 
度 和 稳定 性 。 有 机 质 和 全 和 氮 的 二 阶 微分 建 模 所 得 R? 
分 别 为 0.881 和 “0.956, 为 3 种 建 模 方式 中 最 优 ， 
RMSE 分 别 只 有 2.706 和 0.04$， 能 达到 很 好 的 预测 
效果 。 同 时 ,一 阶 微分 建 模 形式 的 R* 也 能 达到 0.814 
和 0.877， 也 能 得 到 较 好 的 预测 效果 ,而 原 光 谱 建 模 
的 预测 效果 较 差 ， 两 种 土壤 养分 的 预测 模型 及 分 别 
只 有 0.490 和 0.673。 基 于 PLSR 模型 的 光谱 一 阶 微 
分 建 模 和 二 阶 微分 建 模 都 能 很 好 地 预测 土壤 养分 
值 。3 种 建 模 方 式 得 到 的 全 所 的 精度 都 优 于 有 机 质 ， 
说 明 高 光谱 遥感 对 全 所 的 反 演 效果 强 于 有 机 质 。 
2.4 ”人 工 神 经 网 络 建 模 

基于 BP 神经 网 络 分 析 建 立 的 ANN 模型 的 土壤 
养分 实测 值 与 预测 值 的 拟 合 分 布 见 图 4。 选 取 建 模 


波段 作为 自 变 量 ， 以 土壤 养分 值 作为 因 变 量 进行 BP 
神经 网 络 分 析 ， 建 立 原始 光谱 、 反 射 率 一 阶 微分 和 
二 阶 微分 变换 形式 的 分 析 模型 。 由 图 4 可 得 二 阶 微 
分 变换 形式 建立 的 模型 最 优 有机质 和 全 所 二 阶 微 
分 建 模 R? 分别 达 0.886 和 0.984， 优 于 一 阶 微 分 建 模 
的 0.873 和 0.968 以 及 原始 光谱 建 模 的 0.486 和 0.516， 
而 反射 率 二 阶 微分 建 模 的 RMSE 也 分 别 只 有 2.614 
和 0.147。 说 明 二 阶 微分 的 ANN 模型 具有 很 好 的 预 
测 能 力 , 一 阶 微分 模型 也 县 有 较 好 的 预测 能 力 ,而 
原始 光谱 模型 的 精度 和 稳定 性 相对 较 差 。 对 比 有 机 
质 和 全 和 氮 的 3 种 光谱 形式 模型 ,得 到 与 PMLR、PLSR 
模型 类 似 的 结果 ,， 仅 有 机 质 原始 光谱 反射 率 的 稳定 
性 强 于 全 氮 ， 其 他 所 有 指标 均 劣 于 全 氮 , 全 和 氮 的 高 
光谱 反 演 效果 优 于 有 机 质 。 


3 结论 与 讨论 

土壤 光谱 反射 率 是 土壤 特征 研究 的 重要 指标 之 
一 ， 是 土壤 遥感 的 物理 基础 ， 依 靠 其 分 辨 率 高 等 优 
点 已 经 得 到 广泛 应 用 。 探 讨 土壤 光谱 与 土壤 养分 之 
间 的 关系 ,使 用 不 同 数学 方法 建立 模型 并 选择 最 优 
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模型 ， 有 助 于 提高 土壤 养分 遥感 定量 反 演 和 预测 的 
精度 。 本 研究 通过 对 研究 区 土壤 样品 的 光谱 数据 进 
行 去 噪 、 重 采样 等 预 处 理 ， 并 对 土壤 高 光谱 反射 率 
进行 一 阶 微分 和 二 阶 微分 形式 的 变换 ， 对 比 多 元 逐 
步 回 归 分 析 、 偏 最 小 二 乘法 回归 分 析 以 及 人 工 神经 
网 络 分 析 3 种 方法 与 土壤 养分 特征 关系 的 研究 ， 得 
到 以 下 结论 。 

1) 在 3 种 模型 中 , 反射 率 一 阶 微分 和 二 阶 微分 
能 显著 提高 光谱 与 土壤 养分 之 间 的 相关 性 ， 所 得 到 
的 预测 效果 也 均 优 于 原 光谱 反射 率 ， 其 中 二 阶 微分 
的 预测 效果 最 佳 ， 能 很 好 地 预测 研究 区 的 土壤 养分 
状况 。 一 阶 微分 形式 次 之 ,也 能 较 好 地 预测 土壤 有 
机 质 和 全 和 氮 含 量 。 说 明光 谱 变 换 形 式 能 消除 实验 品 
声 等 误差 影响 ， 提 高 光谱 与 土壤 养分 间 的 敏感 性 。 
这 与 前 人 研究 结果 一 致 ， 尤其 是 光谱 的 二 阶 微分 变 
换 中 和 一 阶 微分 变换 中 在 提升 光谱 和 土壤 组 分 之 间 
敏感 性 上 效果 明显 。 

2)3 种 高 光谱 模型 都 能 较 好 地 预测 艾 比 湖 流域 
不 同 土壤 覆盖 类 型 的 养分 状况 。 但 从 建 模 稳定 性 和 
精度 上 而 言 , ANN 模型 最 优 , PLSR 模型 次 之 , SMLR 
模型 最 差 。 其 中 全 和 氮 基 于 ANN 模型 的 二 阶 微分 模型 
效果 最 优 ， 其 决定 系数 达 0.984， 而 均 方 根 误差 只 有 
0.147。 说 明 ANN 模型 在 研究 区 养分 预测 上 具有 很 
好 的 稳定 性 和 高 预测 精度 。 土 壤 光 谱 反 演 是 一 个 复 
杂 的 线性 和 非 线性 相 结合 的 问题 ， 杨 扬 等 上 通过 对 
三 江 源 地 区 有 机 质 的 高 光谱 研究 也 表明 ，ANN 模型 
更 为 适合 土壤 组 分 监测 和 土壤 质量 的 监测 , 但 ANN 
模型 是 否 更 为 适合 其 他 区 域 、 其 他 土质 、 其 他 土壤 
组 分 还 有 待 进一步 研究 。 

3) 在 3 种 模型 的 所 有 光谱 变换 形式 中 , 全 氮 的 
拟 合 效果 均 优 于 有 机 质 。 同 一 模型 中 , 全 氮 预 测 的 
决定 系数 更 高 、 总 均 方 根 误差 更 小 ， 说 明光 谱 对 全 
氮 的 敏感 性 优 于 有 机 质 ， 高 光谱 对 全 氮 的 反 演 效果 
更 好 。 

4) 土 壤 有 机 质 和 全 和 氮 之 间 相 关 性 较 好 ， 达 0.85。 
同时 两 种 养分 与 光谱 数据 的 相关 性 也 很 相似 , 很 多 
反 演 波段 都 很 接近 其 至 一 致 ， 因 此 不 同 土地 覆盖 类 
型 土壤 养分 的 光谱 反 演 能 否 协 同 预测 或 采用 相同 波 
入 进行 反 演 ， 还 有 待 于 进一步 研究 。 


致谢 ”新疆 林业 科学 院 宁 虎 林 、 罗 青 红 、 邹 杰 3 位 
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